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Abstract. The distances to individual wind-driven bubbles such as Planetary Nebulae (PNe) can be determined using expansion 
parallaxes: the angular expansion velocity in the sky is compared to the radial velocity of gas measured spectroscopically. Since 
the one is a pattern velocity, and the other a matter velocity, these are not necessarily the same. Using the jump conditions for 
both shocks and ionization fronts, I show that for typical PNe the pattern velocity is 20 to 30 % larger than the material velocity, 
and the derived distances are therefore typically 20 to 30% too low. I present some corrected distances and suggest approaches 
to be used when deriving distances using expansion parallaxes. 
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1. Introduction 

Distances to individual nebular objects, such as Planetary 
Nebulae (PNe) are important to quantitatively understand their 
structure and evolution, but notoriously difficult to determine. 
Typically uncertainties of a factor of two can be expected from 
so -called s tatistical methods, a review of which can be found 
in Terzian ( 1993). 

Probably the best method for measuring the distance to 
an individual nebula is the 'expansion parallax' method. Here 
the nebular expansion in the sky as measured from images 
taken at different epochs is compared to the radial veloc- 
ity as measured spectroscopically. This method became fea- 
sible with the advent of high resolution imaging, in the late 



1986, 


1989a b; Haiian et alJl993lll995: 


Kawamura & Masson 


1996; 


Haiian & Terzian 1996; Christianto & Seaquist 


1998) 



and nowadays a lso in the optical with the Hubble Space 
Teles cope (HST, iReed et aljfl999t IPalen et alJEool Li et alJ 
2002). 



Clearly the simplest application of this method suffers from 
a number of drawbacks. One complication is the choice for the 
spectroscopic velocity, where different ions often give differ- 
ent values. The other is the assumption of spherical expansion. 
However, authors have been making corrections for the shape 
and aspherical expansion of the measured nebulae, and apply- 
ing such sophisticated templates clearly has made the method 
more useful, in some cases reaching claimed errors as low as 
10-20%. 
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Still, the application requires a reasonably close or partic- 
ularly rapidly expanding nebula and to date has only been ap- 
plied to a limited number of PNe. Table 1 list all cases with 
a well determined expansion parallax distance, leaving out the 
cases where only upper limits were found. 

Among the corrections needed to accurately use the expan- 
sion parallax method there is one which to date has received lit- 
tle attention in the literature, even though it is quite basic. This 
is the fact that the expansion velocity as measured in the sky 
is a pattern velocity, whereas the spectroscopically measured 
velocity is a material velocity, and the t wo are generally not 
the s ame. This effect wa s touched upon by Mart en et alJ (fl993) 
and lSteffen etaP il997h . but has not been taken into account in 
any of the published expansion parallax distances. 

The extreme case would be that of an R-type ionization 
front making its way through a stationary medium, and clearly 
here the expansion parallax method becomes useless, as has 
been mentioned by various authors. Luckily, this situation is 
thought to be rare in PNe. However, also in the more com- 
mon cases of shock fronts or slower moving D-type ionization 
fronts, the two velocities will differ. Although typically not by 
much, the effect is systematic, not random, and it should there- 
fore be taken into account when using the expansion parallax 
method. Especially as the measurements are becoming more 
accurate over time, this factor can no longer be neglected. 

In this paper I calculate the magnitude of the discrepancy 
between the two velocities for shocks (Sect. 2), ionization front 
structures (Sect. 3) and so-called 'shells' (Sect. 4). In Sect. 5 
I apply the derived corrections to the published results of the 
expansion parallax method and suggest strategies to optimize 
the expansion parallax method. 
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Table 1 - Expansion parallaxes 



PN 


^spectro (km S ) 


Distance (kpc) 


R 


BD+30 3639 


22 ± 4 


2.8+ 12 


r 


BD+30-3639 


22.0 ± 1.5 


2.68 ± 0.81 


r 


BD+30 3639 


22 ± 4 


1.5 ± 0.4 


r 


BD+30 3639 


25.6 


1.2 ± 0.12 





IC 2448 


17.9 ± 0.3 


1.38 ± 0.4 





NGC 3242 


26 ± 4 


0.42 ± 0.16 


r 


NGC 6210 


23 ± 5 


1.57 ± 0.40 


r 


NGC 6543 


16.4 ± 0.16 


1.00 ± 0.27 





NGC 6572 


14 ± 4 


1.49 ± 0.62 


r 


NGC 6572 


14 ± 4 


1.2 ± 0.4 


r 


NGC 6578 


19.2 ± 0.5 


2.00 ± 0.5 





NGC 6884 


16.6 ± 0.4 


2.20 ± 0.8 





NGC 7027 


21 


0.94 ±0.2 


r 


NGC 7027 


17.5 ± 1.5 


0.88 ±0.15 


r 


NGC 7027 


17.5 ±1.5 


0.703 ± 0.095 


r 


NGC 7662 


21 ±7 


0.79 ± 0.75 


r 


VY 2-2 


19.5 ±0.4 


3.6 ±0.4 


r 



Radio/Optical Reference 



Masson (1989b) 
Haiianetal. (1993) 
Kawamura & Masson (1996) 
Li et al. (2002) 
Palen et al. (2002) 
Haiian et al. (1995) 
Haiianetal. (1995) 
Reed et al. (1999) 
Haiianetal. (1995) 
Kawamura & Masson (1996) 
Palen et al. (2002) 
Palen et al. (2002) 
Masson (1986) 
Masson (1989a) 
Haiianetal. (1993) 
Haii an & Terzian (1996) 
Christianto & Seaauist (1998) 
* This is the spectroscopic velocity used by the authors, which can correspond to a shock or ionization front velocity 
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Fig. 1. Sketch of the shock configuration and the definition of 
the various velocities. 



2. Shock waves 

Since PNe are basically wind-driven bubbles, they contain 
shocks. Numerical hydrodynamic modelling has shown that 
in round or elliptical nebulae the inner bright rim is nor- 



mally associated with the shock wave being driven into the 
surrounding material by the stellar wind. Key papers de- 
scribing the radiation-hydrodynamic evolution of P Ne are 
iMarten & Schoenbernerl(ll99ll):lMellemal(ll994lll995l) . and a 
review is presented in Schonb erner & St effen (2003). The ap 



plication to an individual PN is shown in lco^adTetaD J200C 
Whenever I mention numerical simulations in what follows, I 
refer to these papers. 

The jumps of density and velocity across a shock are given 
by the Rankine-Hugoniot conditions. These are usually given 
for the reference frame of the shock, in which case the velocity 
jump is 



{ 1 + l)M 2 
( 1 -\)M 2 + 2 



(1) 



which for infinite Mach number M. and an adiabatic index 7 = 
5/3 gives the classical value of 4. The Mach number of the 
shock is given by 



M 



a 



(2) 



where ciq is the sound speed in the pre-shock gas, given by 



a 



hkT 



(3) 



with k the Boltzmann constant, To the gas temperature, fi the 
mean molecular weight, and toh the mass of atomic hydrogen. 

However, in the stellar frame, the shock has a velocity 
which we will call u s and the pre- and post-shock velocities 
are given by w .i = v o,i ± u s, see Fig-E 

The expansion parallax method measures u s from the an- 
gular expansion of the shock front, 9 = u s /D, and u\ from the 
spectroscopy, and derives the distance D from the ratio of the 
two, assuming u s = u\. From the expression above it is imme- 
diately clear that u s and u± are not equal. If the PN woud be 
expanding into complete vacuum the two velocities would be 
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identical, but in the case where it is moving into a surrounding 
medium, it is more proper to view the expansion as a distur- 
bance or wave travelling through the medium, raising its den- 
sity and velocity. The speed at which the wave travels is not the 
same as the velocity to which the material is accelerated by the 
passage of the wave. 

Let us consider the ratio 1Z of the two velocities, u s /u±. 
The distance to the PN can then be expressed as 1Zu\ / 6. Using 
the relations between Uo,i and Vo,i together with Eqs.n an d|2] 
(realizing that in our shock frame both vq and v% are negative) 
gives 



K 



(7 + l)Mu a + (7 + l)M 2 a 
(7 + l)Mu + 2{M 2 - l)a 



(4) 



The limits of this ratio are (uq + ao)/uo for M. — > 1 and 
(7 + l)/2 for M. — > 00. This shows that only for isothermal 
(7 = 1) hypersonic shocks the ratio tends to one for high 
Mach numbers. Choosing values for the pre-shock velocity uo, 
pre-shock sound speed ao, and shock Mach number M. gives a 
value for 1Z. The material velocity u\ can then be found from 



Ui = Uq + 2a 



M 2 - 1 



(5) 



Reversely, given values for uq, U\, and ao, the shock's Mach 
number can be found from 



M = (7+ i)(m - up) + Vh + i) 2 K - m) 2 + 16^ 

4a 

which can then be used to derive 1Z using Eq.|4] 

In Fig. |2 1 plot the ratio 1Z as function of u\, for six pre- 
shock velocities u in the range 1 to 25 km s _1 , using a 7 
value of 5/3 (monatomic ideal gas), and a = 15 km s _1 , 
the adiabatic sound speed for an ionized gas at 10 4 K. Most 
PNe have high densities and relatively slow shocks, so that the 
shocks are expected to be isothermal. Using 7 = 1 is therefore 
appropriate, and in Fig. [3] I show 1Z as function of u\ for this 
isothermal case, where I have used an isothermal sound speed 
of 1 1.7 km s _1 , valid for an ionized gas of cosmic abundances 
at a temperature of 10 4 K. As expected from the limits derived 
above, the ratio is largest for low velocities. I should note that 
the ratio does depend on the choice for the sound speed a . A 
value corresponding to an electron temperature of 10 4 K is typ- 
ical for PNe, but the temperature can be both higher and lower 
than this, ranging from 5000 to 15,000 K. Lower temperatures 
give lower values of 1Z. 

Looking at the data in Table 1, it is clear that for these 
PNe, we are never in the very high Mach number regime. Using 
Fig-EJone can see that for the observed ranges of velocities and 
for 7 = 1, the typical ratio 1Z is between 1.3 and 1.5, although 
it can be as high as 1.8. Obviously this falls outside the formal 
errors of the method and becomes as important an effec t as th e 
geometric corrections applied for example bv lLi et alJ (|2002). 
Furthermore, the factor is always larger than 1, so it does not 
make sense to add a 20-30% extra uncertainty to the distances. 
Rather, the distances should be scaled up by 30% and an extra 
error of ~ 10% added to it. Figure[5]can be used to estimate 
the magnitude of this effect for individual PNe, which I will do 
in Sect. 5. 




10 20 30 40 

Uj (km s" ) 



50 



Fig. 2. The ratio TZ(— w s /iti) as function of u% for the case 7 = 
5/3 (no cooling), and (adiabatic) sound speed ao = 15 km s _1 . 
The six curves correspond to different values for the pre-shock 
velocity uo: 1 (highest curve at ui = 50 km s -1 ), 5, 10, 15, 20, 
and 25 (lowest curve) km s _1 . 




10 20 30 40 

Uj (km s" 1 ) 



50 



Fig. 3. The ratio 7Z (= u s /u\) as function of u\ for the 
isothermal case (7 = 1), and (isothermal) sound speed ao = 
11.7 km s _1 . The six curves correspond to different values for 
the pre-shock velocity uq: 1 (lowest curve), 5, 10, 15, 20, and 
25 (highest curve) km s _1 . 

The conclusion then is that if the expansion parallaxes are 
measured from the shocked component of the PN, the distances 
should be multiplied by a factor 1Z given by Eq. |4j where un- 
certainties in the velocity of the material into which the shock 
is expanding adds a ~10% error. I will come back to this point 
in Sect. 5. 

3. Ionization fronts 

Ionization fronts are either R- or D-type. The basic difference 
is the speed of the front which for R-type fronts is higher than 
twice the sound speed in the neutral gas, giving the gas no time 
to react to the presence of the ionization front, wheras for D- 
type it is lower than twice the sound speed in the neutral gas. 
The result is that D-type fronts actually consist of a combina- 
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10 12 14 16 18 20 

u 2 (km s" 1 ) 

Fig. 4. The ratio 1Z (= Ui/u^) as function of u-i for a weak 
R-type ionization front running into a 10 km s _1 wind. The 
weaker the front the closer the ratio gets to one. 




u 2 (km s" 1 ) 



Fig. 5. The ratio 1Z (— Ui/112) as function of U2 for a weak D- 
type ionization front running into a wind with velocity u\ of 
(from left to right) 20, 25, and 30 km s _1 . 



tion of an ionization front and a preceding shock front. A good 
overvi ew of basic ionization front theory can be found in Shu 
Jl992h . 

The numerical simulations for the formation of PNe which 
studied the effects of ionization fronts in detail, showed that 
although R-type fronts occur, they never persist for long, and it 
would be highly unlikely to observe one in action. When they 
occur, they are of the so-called weak type, i.e. supersonic with 
respect to both the neutral and ionized flow regions. 

Weak D-type fronts (moving subsonically with respect 
to both the ionized and neutral regions) are more commonly 
found in the models, and persist for a longer time, and could be 
present in a number of observed PNe, although it is generally 
hard to prove this. The simulations also show that the presence 
of a D-type front sets up a disturbance in the gas which persists 
long after the front itself has disappeared, creating the shells 
which I discuss in the next section. 

Basic ionization front theory shows that the velocities v\ 
and V2 on either side of the front are related by 



- = A (4+vf ±[(al + vir-4alvl] 

V2 Ztto ^ 



1/2 



(7) 



The indices 1 and 2 refer to the neutral and ionized sides, re- 
spectively. Just as in the previous section, v refers to velocities 
in the frame of the discontinuity, and u to velocities in the stel- 
lar frame. If the front moves with a velocity it; in the stellar 
frame, the velocities on either side of it are U1.2 = t»i,2 + Ui. 
For the weak R-type fronts one has to use the — sign in Eq. 
and for the weak D-type fronts the +. The velocity v\ has to 
be larger than vb, = a 2 + (a| ~~ ct i) 1 ^ 2 f° r R-type fronts and 
smaller than v d = a 2 — {a\ — a]) 1 / 2 for D-type fronts. 

In Fig. |^I show the ratio 1Z, now defined as Ui/u 2 , against 
u 2 for weak R-fronts running into a medium with ui 
1 km s . For the sound speeds I chose a\ = 1 km s _1 
and 0,2 = 10 km s _1 , and tests show that the ratio 1Z is hardly 
sensitive to this choice. As can be expected this ratio is large for 
very fast moving fronts, and approaches 1 for the slower mov- 
ing ones. Because of the large range of ratios, it would seem 



dangerous to apply any type of correction in case of R-type 
ionization fronts. 

In Fig.[5]l show the same ratio 1Z against 112 for weak D- 
fronts, a situation which is more likely to occur in real PNe. 
Normally the shock front which precedes the D-front will have 
accelerated the AGB wind to a higher velocity. The numeri- 
cal simulations show this velocity to be typically around u\ = 
20 km s _1 . 1 plot the ratio for three values of u%. Interestingly, 
around the typically observed velocities, the correction is of the 
same order of magnitude as for the shock waves in Sect. 2. 

The full s olution for a D-type front would eliminate the 
choice for u\. Ishuldl992l) showed how this can be done for 
a constant flux of ionizing photons impinging on a static and 
constant density environment. Attempts to do this for a stel- 
lar wind type environment have been only partly successful. 
Giuliani (1989) found a self-similar solution which requires the 
flux of ionizing photons to be time dependent and fall off as 
which does not apply to most PNe. 

Masson (1986, 1989a b) applied a series of corrections to 
his determination of the distance to NGC 7027, BD+30 3639 
and NGC 6591, which are partly related to the difference be- 
tween the material velocity (which he referred to as the "bulk 
velocity"), and the (pattern) speed of the ionization front. The 
magnitudes of these corrections are actually similar to the ones 
found above, 1.1 to 1.2. Perhaps somewhat confusingly he ap- 
plies these corrections to the measured angular expansion rate, 
so that the final figures he quotes for 9 are not actually the ones 
measured. He derives this correction from the fact that as the 
ionized shell expands, the density will go down, reducing the 
number of photons used up in recombinations, and allowing the 
ionization front to expand. This is a different approach to the 
one used above and does not use the jump conditions. However, 
the basic idea is still that the ionization front expands faster 
than the material flow of the gas, so the corrections are related. 
In not considering the jump conditions he implicitly assumes 
a weak R-type ionization front with a small density jump, and 
the correction corresponds therefore to the one from Fig.|4] The 
numerical simulations as well as the observed velocities show 
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f 1 .30 




Uj (km s" 1 ) 



Fig. 6. The rat io 1Z (= Ms/mi) as f unction of i t i for t he shell so- 
lutions from lChevalierl 119971) and lShu et all J2002h . The solid 
line shows the solutions for a photo-ionized wind, the dotted 
line the same but with the influence of a wind included. The 
dashed line shows the result for photo-ionized winds with dif- 
ferent density laws steeper than r~ 2 . 

that this is not correct, although the correction factors come out 
at similar values. Note that Masson is the only author to actu- 
ally apply any type of correction. In Sect. 5 I will comment 
some more on this. 

The conclusion from this section is that it is more difficult 
to find a general solution for ionization fronts, which makes 
them less suitable than shock fronts for applying the expansion 
parallax method. However, the results in this section allow for 
corrections to be made, and below I will illustrate how these 
corrections could be applied. 

4. Shells 

Numerical simulations show that the density disturbance in- 
duced by the D-type ionization front will persist and can be 
identified with the observed attached shells aro und the br ight 
core nebula, such as for example in NGC 3242. ICorradi et alJ 
(2003) compiled a list of PNe with haloes and shells, which can 
serve as a reference. We follow their nomenclature and refer to 
these structures as shells. 

Analytica l models f or the expa nsion of such she lls are 
available in IChevalierl J 19971 and IShu et alJ (E)02), even 
though the latter had a different application in mind. Using 
their results, it is possible to extract a similar ratio 1Z of pattern 
speed over material speed, where the pattern speed is now the 
movement of the edge of the shell. Since these models actually 
use the isothermal shock conditions, the results are basically 
identical to ones from Sect. 2. However, I lis t them separately , 
sinc e shells are com monly observed and the lChevalien (|1997) 
and lShuetalJ 12002) results are the full solutions including the 
acceleration of the AGB wind after ionization. 

For an AGB wind of co nstant mass loss and velocity, the 
results from lChevalierl(ll997l) show that 7Z is between 1.10 and 
1 .34. 1 plot this ratio against u\ in Fig. [6] together with the ratio 
from his solution including the effect of a stellar wind, which 
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Table 2 - New distances 



rlN 


Distance (kpc) 


BD+30 3639 


1.3 ±0.2 


IC 2448 


2.07 ±0.62 


NGC 3242 


0.55 ±0.23 


NGC 6543 


1.55 ±0.44 


NGC 6578 


2.90 ± 0.78 


NGC 6884 


3.30 ± 1.24 


NGC 7027 


0.68 ±0.17 


NGC 7662 


1.19 ± 1.15 


VY 2-2 


4.68 ± 1.20 



gives somewhat lower values. The ratio for an AGB wind den- 
sity falling off steeper (p oc r~ a with a > 2. IShu etalJl2002l) 
are also shown. 

The numerical simulations also show that when an attached 
shell forms, its velocity structure is that of a rarefaction wave, 
with a positive outward gradient. This means that the bright rim 
expanding into this will be moving into an area with a velocity 
lower than the original AGB wind velocity. The implication for 
the expansion parallax method is that when using the rims of 
PNe with attached shells a value of uq lower than 10 km s _1 is 
more appropriate when determining the correction factor from 
Sect. 2. 

The conclusion from this section is that for attached shells 
a correction of around 20% is needed. 

5. Implications for distances and the method 

The results presented thus far show that for the published dis- 
tances, the results should typically be scaled up with a factor 
1.3 ± 0.2, which interestingly enough eliminates the system- 
atic discrepancy between the expansion paralla x distances and 
statistical distances noted bv lPalen et ail (2002). In this section 
I will go through the list of PNe for which distances have been 
measured, consider corrections for each case, and give some 
suggestions on how to improve the usage of the expansion par- 
allax method further. For calculating the corrections I use the 
values 7=1 and gq = 11.7 km s _1 for the cases of shocks 
(corresponding to Fig.|5|l, a\ = 1 and ai = 10 km s _1 for 
ionization fronts (corresponding to Fig.|5Jl. The first part of this 
section is intended not only to derive new distances, but also to 
illustrate how to derive and use the correction factors. 

5. 1. New distance estimates 

Going through the list there are a number of PNe which con- 
tain bright rims and attached shells, and the rim can therefore 
be considered to be bounded by a shock. These are NGC 3242, 
NGC 6578, NGC 6884, NGC 7662, and IC 2448. Making con- 
servative assumptions about the value of uq (13 ± 12 km s _1 ) 
and using the results from Sect. 2, the distances to these PNe 
should be increased by factors 1.3±0.2, 1.45±0.15, 1.5±0.15, 
1.5 ± 0.3, and 1.5 ± 0.1 respectively. The uncertainties are due 
partly to the uncertainties in the reported spectroscopic veloci- 
ties and partly to the uncertainty in the value of uo ■ The latter 
could in principle be reduced if one assumes that the velocities 
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in the shells follow the numerical models, in which case uq will 
be low (< 10 km s _1 ). I list the new (conservative) values in 
Table 2. 

NGC 6543 is a more complex nebula. The Chandra results 
show that it contains a wind-driven bubble, but the ionized 
material surrounding this bubble does not resemble a standard 
shell. This complicates the choice for Urj, but within reasonable 
limits the correction factor is 1.55 ± 0.15. 

Three PNe in Table 1 are candidates for the presence of an 
ionization front: BD+30 3639, NGC 7027, and Vy2-2. 

Of these BD+30 3639 is a single shell, and relatively young 
PN with molecular emission around it, although the deeper 
exposures also seem to show an ioni z ed component aro und 
the bri ght shell CS ahai & TraugeJl998l) . lKawamura & Massonl 
i 19961) used a detailed analysis of the expansion of this object 
in order to derive the distance, assuming we are observing an 
ionization front. Part of this analysis is a correction factor of 
1 . 19 to account for the difference between the material and pat- 
tern velocities, as mentioned in Sect. 3. Using the results from 
Sect. 3 for a weak D-type front, I estimate a correction factor of 
1.2 ± 0.2 assuming a velocity of 25 ± 5 km s~ x for the shocked 
neutral g as and a spectroscopic veloc ity of 22 ± 4 km s _1 (as 
used by iKawamura & Massonlll9 96). confirming their result 
for the distance. 

Li et al. used optical HST data, better spectroscopic data 
and a correction for the ellipticity of the nebula to arrive at 
a smaller distance. Using their preferred value for the spec- 
troscopic velocity, 25.6 km s _1 , I find a correction factor of 
1.1 ± 0.1 if we are seeing an ionization front running into a 
shocked component with a velocity of 25.6 - 30 km s _1 , and 
a correction of 1.3 ± 0.1 if we are seeing a shock. It therefore 
seems reasonable to increase the distance by 10% and to dou- 
ble the uncertainty on that numb er: 1.3 ± 0.2 kpc, b ringing it 
closer to the value from lKawamura & Masson| fl996) 

NGC 7027 is a more complex PN where an ionization 
front has to be present, given the large amounts of dust and 
molecules surrounding this object. The elliptical shell seen at 
radio wavelengths (and partly obscured at optical wavelengths) 
could be the io nization front and is analyzed as such by M assonl 
(1986, 1989a). To correct for the discrepancy between spec- 
troscopic and angular velocities, he uses a factor of 1 .2 in the 
first and 1.15 in the second paper (note that the second paper 
quotes 1/1.15=0.83 as the correction factor). Further correc- 
tions for the decrease of the radio flux in the second paper ac- 
tuall y largely cancel th e effect of this factor. 

iBains et aD d2003l) analyzed high resolution optical long 
slit spectra of NGC 7027 and from the [O III] line derived an 
equ atorial expa nsion velocity considerably lower than use d by 
lMassonldl989al) . namely 13 ± 1 km s' MCox et ail J2002I) re- 
port K-band imaging and spectroscopy for NGC 7027. Their 
best fitting model has an equatorial expansion velocity in Br7 
of about 13 km s _1 (Cox & Huggins, private communica- 
tion; this number is not given in the paper), consistent with 
the [O III] value. For the molecular emission (H2) the same 
authors find an equatorial velocity of ^15 km s -1 . Allowing 
all possible values for the velocity of the neutral material u\ 
(13 - 23 km s _1 ), the correction factor would be 1.4 ± 0.4. 
Assuming that the H2 expansion is indicative of the value of u\ 



(13 - 17 km s ), the ratio becomes 1.2 ± 0.2. If the ionized 
shell is actually bounded by a shock, the results of Sect. 2 show 
that the correction factor would be 1.75 ± 0.15 for all allowed 
values of the pre-shock velocity u§. 

Tak ing an a ngular expansion of 4.84 ± 0.82 mas yr _1 
(Masson 1989a[ this includes the correction for flux variations 
0.47 ± 0.47), but not the correction factor of 0.83), combin- 
ing this with an equatorial velocity of 13 km s , and apply- 
ing a correction of 1.2 ± 0.2, I arrive at a new distance of 
680 ± 170 pc. Note that this error is rather optimistic in view 
of the wide range of correction factors mentioned above. 

Vy2-2 is a compact PN which is hard to categorize. Also 
here one may have to be aware of ionization fronts since this is 
a fairly young, low excitation PN. Assuming either an ioniza- 
tion front or an isothermal shock, the reported distance should 
be scaled up by 1.3 ± 0.3 

Two PNe in the list are hard to categorize in the standard 
wind-blown bubble plus photo-ionization scheme: NGC 6210 
and NGC 6572. 1 therefore will not suggest any corrections for 
these, although shocks could very well be present. For NGC 
6210 this would imply a 30% increase in the reported distance, 
but for NGC 6572 the correction could become very large ( 1 .4- 
1.8) due to the low value of the reported spectroscopic velocity. 

5.2. Improving the method 

The correction factors derived in this paper depend on a num- 
ber of parameters. As already indicated, due to the isothermal 
nature of the slow shocks in PNe, the preferred value for the 
adiabatic index 7 is 1 . The sound speed in the ionized medium 
(ao) I have taken to be 11.7 km s _1 in Fig. [5] which is the value 
for an electron temperature of 10 4 K. If the PN is known to 
have a particularly high or low electron temperature the figure 
should be recalculated for the appropriate sound speed, since 
the values for 7Z unfortunately do depend on the choice of the 
sound speed. The appropriate sound speed to use here is the 
isothermal sound speed. 

Figure [3] or its equivalent can then be used to find the 
correction factor for the range of values of uo and u\ which 
seem reasonable. Observational data on uq is scarce, but if it 
cannot be further constrained, 10 ± 10 km s _1 should cover 
most cases. For u\ the question is which spectroscopic ve- 
locity to pick. Ultimately, for an individual PN this question 
can only be answered through detailed (photo-ionization or 
hydrodynamic) modelling of the PN, assuming it is actually 
possible to produce a unique model for it. The approach of 
iGesicki et alJ Jl998t) . who derive spatial velocity profiles using 
photo-ionization models, may be useful here. Generic hydro- 
dynamic modelling may help establish what would be a good 
choice in general cases and I understand that a project to do 
this is under way (Schonberner, private communication). Until 
that time the choice for u\ will introduce uncertainties in the 
method, as it has always done. It is important to realize that the 
jump conditions are valid just before and after the front, so the 
best choice would be a velocity as close as possible to the front. 

Round or elliptical PNe with attached shells are particu- 
larly suitable for the expansion parallax method. Their inner 
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bright rim is bounded by a shock, and the presence of a sur- 
rounding ionized shell shows that this is not an ionization front. 
For these almost round PNe, long slit spectroscopy can give a 
good indication of the spatial and velocity structure, since it 
will be quite similar in the directions perpendicular and par- 
allel to the line of sight. This way a distinction can be made 
between the velocity of the rim and the shell (the latter can be 
both higher and lower than that of the rim). If possible it is 
best to choose species which lie closest to the outer edge of the 
rim so as to avoid any velocity variations within the rim. For 
these rims numerical models show that they are expanding into 
a region of low velocity, typically lower than that of the orig- 
inal AGB wind. Choosing a value for uo of around 5 km s -1 
is a reasonable guess. In any case, the corrections as derived in 
Sect. 2 only weakly depend on this choice. Using the measured 
spectroscopic velocity, the ratio 1Z can be determined. 

It would be interesting to derive an independent distance 
using the attached shells, if available. If the long slit spec- 
troscopy provides a velocity for the shell, then the corrections 
from Sect. 4 can be used to find a distance. Models show that 
the velocities of the shells are not as strongly position depen- 
dent as those of the inner rims, something which is reflected 
in the only mild ellipticities found observationally. This may 
make the outer edges of shells actually better suited for dis- 
tance determinations with the expansion parallax method. 

It is best to avoid areas and features well away from where 
a shock is suspected. At least for a shock we can make a con- 
nection between the pattern and material velocities. Away from 
these discontinuities, the two will also differ, but there is no 
way to know how. 

Bipolar PNe display a larger range of expansion velocities, 
which means that it is harder to correct for inclination effects. 

II j i 
Palen et al. ( 2002) state that the method is unusable for extreme 

bipolars, which seems too pessimistic. However, a thorough 
understanding of the dynamics of the complete PN is essen- 
tial to be able to apply the method to these PNe, and in most 
cases this information is not available. It would also be valu- 
able to have more radiation-hydrodynamic modelling of such 
systems. 

6. Conclusions 

Pattern velocities and the material (or bulk) velocities in a gas 
are not necessarily the same. For discontinuities, such as shocks 
and ionization fronts, the relation between the two can easily 
be derived using the jump conditions across them. The pattern 
velocity is then found to be always higher than the material 
velocity. Since measuring the expansion of PNe in the sky is 
mostly done using sharp edges, which are associated with ei- 
ther shocks or ionization fronts, a correction should be applied 
before calculating the distance from the ratio of the two ve- 
locities. This correction is typically larger for velocities of the 
order one or two times the sound speed in the ionized material, 
which is actually what is measured in most PNe. 

Not using this correction will systematically underestimate 
the distances to PNe. For the sample of PNe to which the ex- 
pansion parallax method has been applied successfully, the cor- 
rection factors are around 1.2 to 1.3 for both shocks and ioniza- 



tion fronts. Applying the corrections given in this paper should 
lead to improved distance determinations to PNe. 
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